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ABSTRACT 

The wave function of the universe is evaluated by using the Euchdean path 
integral approach. As is well known, the real Euclidean path integral diverges 
because the Einstein-Hilbert action is not positive definite. In order to obtain a 
finite wave function, we propose a new regularization method and calculate the 
wave function of the Friedmann- Robertson- Walker type minisuperspace model. 
We then consider a homogeneous but anisotropic type minisuperspace model, which 
is known as the Bianch type I model. The physical meaning of the wave function 
by this new regularization method is also examined. 
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Quantum cosmology is one of the most fascinated subjects in modern physics. 
In quantum theory, the wave function characterizes each theory. Therefore the 
main purpose in quantum cosmology is to evaluate the wave function of the uni- 
verse. As is well known, the wave function of the universe is obtained as the solution 
of the Wheeler-DeWitt second-order functional differential equation. This equa- 
tion, however, is very difficult to solve even in minisuperspace models, although 
in minisuperspace models the Wheeler-DeWitt equation is reduced to an ordinary 
differential equation. In addition, we have to know the initial condition of the 
universe which is hard to understand completely. Another approach, which gives 
the wave function by the Euclidean path integral, has been proposed by Hartle 
and Hawking ' . They claimed that the path integral should be performed over 
regular and closed four- dimensional manifolds. This idea is called the boundary 
condition of no-boundary. The four-dimensional manifold must be Euclidean to 
avoid the cone singularity at an initial point. Their proposal solves the initial con- 
dition problem but the real Euclidean path integral does not converge because the 
Einstein-Hilbert action is not positive definite. The question is therefore whether 
we can construct the Euclidean path integral approach which is free from a diver- 
gence problem. 

A first way to avoid this difficulty was proposed by Gibbons, Hawking and 
Perry ' . They suggested that the path integral could converge if we performed 
the conformal rotation. This proposal, however, has several objections . Hartle 
proposed another path integral approach. He suggested that the path integral 
should be taken along the steepest-descent path in the space of complex four- 
metrics. Following the Hartle's alternative proposal, we can make the Euclidean 
path integral convergent by changing an integration contour from a real axis to 
a curve on a complex plane. Halliwell and Louko calculated the wave function of 
the universe following this method ' ' . The integral along this complex curve 
is so difficult to perform that they evaluated the wave function in a saddle point 
approximation . Moreover the physical interpretation of this wave function is 
hard because the path integral is defined on the complex four-manifold 



The contour changing method seems not to complete and it is worthy investi- 
gating a possibihty of another regularization method. 

In this paper, we propose a new regularization method and evaluate the wave 
function in minisuperspace models. We first propose a new regularization method 
and calculate the wave function of the Friedmann- Robertson- Walker type minisu- 
perspace model. We then consider a more general case, i.e. a homogeneous but 
anisotropic type minisuperspace model, which is known as the Bianch type I model. 
We give a physical interpretation of the wave function by this new reguralization 
method. 

We first calculate the wave function of the closed Friedmann- Robertson- Walker 
type minisuperspace model ' . The metric of the Friedmann-Robertson- Walker 
type is 



ds^ = a^l^'-L^dr^ + q{r)dnl}, 



where a^ = 2G'/37r and dVi^ is the metric on the unit three-sphere. The Einstein- 
Hilbert action with the rescaled cosmological constant A then becomes 



1 r . ..r g2 



/bW] = 2|, rfriV[-^ + Ag-lJ. (2) 

We can rewrite this action to the Hamiltonian form as follows. 

I\p{T),q{T)] = J rfr|M-^(-V_Ag + l)}, (3) 

where T = N ■ {t" — t') and p = —q/4N. The Wheeler-DeWitt equation can 
be easily obtained by replacing p in the Hamiltonian with —d/dq, therefore the 
Wheeler-DeWitt equation in the Friedmann-Robertson- Walker type becomes 



i7^ = - 

2L 



rf2 



^ = 0. (4) 
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The wave function of the universe can be obtained if one solves this equation . As 



is stated before, we can estimate the wave function by means of the path integral 
^j = I dT f DpDqe^^^P^'^^ = I dT^{T), (5) 
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instead of solving the Wheeler-DeWitt equation . The Wheeler-DeWitt equation 
will be a guide for the path integral to define its contour. As is easily seen, the 
integration with respect to T is not easy to evaluate. We first consider the p and 
q integration . To perform the p and q integration, we divide the interval T into 



L + 1 pieces. Then ip{T) is 

L L 



dpi 

2n 



(6) 



^(T)= hm [fldqifl 

where AT = T/(L + 1), Ti = i ■ AT, qi = q(Ti), pi = p(Ti). The p integration is 
carried out immediately because of the Gaussian form with respect to p. In order 
to perform the q integration, we use a following formula which is obtained by the 
mathematical induction 

= (n + 1)1/2 U"'^ (7) 

^^Y>{^^{h- af ~ af3['^a+'^h- ^{n + l){n + 2)(3\]. 

We can evaluate the p, q integration completely and the result is 

where 

V ; 24 1 4 2^ 8T ^ ' 

In Ref.[8], Halliwell and Louko divided the variables p{T) and q{T) into the classical 



part and the quantum one. On the other hand, we can calculate the p, q integration 
directly by using above induction (7). 

We now evaluate the T integration. T is the rescaled lapse function and the 
range of the T integration must preserve the time reparametrization invariance. 
The T integration should be performed from — oo to +00 naively. This integra- 
tion, however, diverges if we take the contour along the real axis. From this reason, 
a regularization is necessary to remove this divergence. In Ref.[8, 9, 10], the inte- 
gration is regularized by means of changing a contour of the T integration and, by 
using a saddle point approximation, they evaluated the wave function. While the 
contour changing method is a very natural regularization, the full integration is 
very difficult. It seems to be hard to study the wave function beyond a saddle point 
approximation. We adopt another strategy, i.e. instead of changing the contour 
directly, we extend the integrated function ip{T) in (8) such that the T integration 
will converge. 

We introduce a complex parameter a and define the wave function as an ana- 
lytic function with respect to a as follows (up to a normalization factor). 

where the rescaled lapse function T remains real and the complex parameter a must 
be determined in such a way that the T integration will converge. The strategy is as 
follows. We temporally fix a as a constant complex number during the integration 
with respect to T. After the integration, we take the analytic continuation a — * 1 
and get the wave function. 

Let us perform the above strategy and obtain the wave function of the Friedmann-| 
Robert son- Walker minisuperspace. In this case, unfortunately, we cannot carry out 
the integration completely. However we can use the Wheeler-DeWitt equation and 
obtain the wave function, which is equivalent to calculating the integration. This 
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method was proposed by Halliwell and Louko . Apparently the wave function is 



symmetric with respect to g' = qir') and q" = q{T"). In addition, the wave function 
satisfies the Wheeler-DeWitt equation at q' and q". By rescahng T -^ A~^/^2^/^T 
and introducing a new variable z = (1 — Ag)/(2A) ' , the wave function as a solution 
of the Wheeler-DeWitt equation will be expressed as Airy function's products, 

^ =aAi{z")Ai{z') + hBi{z")Bi{z') 

(11) 
+ c[Ai{z")Bi{z') + Bi{z")Bi{z')\ . 

In order to get coefficients {a,b,c), we take z' = z" = z and expand (11) with z 
around z = 0. On the other hand, (10) will be 

where T is complex in this representation. We also expand this with z around 
2; = 0, and get coefficients (a, 6, c) by comparing their coefficients of z. We suppose 
that z is very small value and we need not 0(z"^) terms. 

As is stated before, the complex parameter a must be so determined that the 
integration will converge. This statement requires that i?e[(«T) ] is positive for 
\T\ ^ 1, hence there are three a's areas. In each case, we get the wave function. 
The results are 

^ = Ai{z")Ai{z'), (13) 

^ = Ai{z")Bi{z') + Bi{z")Ai{z'), (14) 



^ =Ai{z")Ai{z') - Bi{z")Bi{z') 

± t [At{z")Bt{z') + Bt{z")Bt{z')] , 



(15) 



(16) 



^ =3Ai{z")Ai{z') + Bi{z")Bi{z') 

± i [Ai{z")Bi{z') + Bi{z")Bi{z')] , 

where the overall factor is neglected. The wave function (13) has the same form 
of the Lorentzian one in Ref.[7]. In this case, the a's area is constrained on a pure 



imaginary axis in our method. Therefore, up to an overall factor, the path integral 
form of the (13) 's is the same with the Lorentzian one. We can show that the wave 
function becomes the (13) 's form without using the Wheeler-DeWitt equation, by 
following the inverse of the Halliwell's calculation in Ref. [7] appropriately. This fact 
gives us the insight that the singularities in (8) or (9) is caused for the undesirable 
integration. We can eliminate singuralities if we consider the products of Airy 
functions. This is why we take z' = z" = z. On the contrary to the (13) 's case, 
the path integral form of the (14)'s , (15)'s and (16)'s cases can not be calculated 
directly. Therefore we use above insight and we can calculate the wave function for 
other three cases, those are (14), (15) and (16). In our regularization method, we 
can follow this insight easily because we have only turned up the straight contour 
line by using the complex parameter a. Of course we must take care of a branch cut 
caused by the T"^'^ term in the complex T plane in (12). We consider both cases 
that the contour passes through the cut or not. From this reason, we get six results 
with three a's areas. On the other hand, in the contour changing method, the 
contour is a complicated curve. It is very difficult to perform this insight beyond a 
saddle point approximation and to distinguish whether the contour passes through 
the cut or not. Therefore the results in Ref. [8] are the part of our results. 

Of course we can take another strategy to remove a divergence problem. If we 
extend the integrated function as 

where j3 is some constant, then \E' also converges. In this representation, we also fix 
a as Re[{aT)'^] is positive for \T\ ^ 1. After the integration, we take the analytic 
continuation a ^ 1. We also calculate the wave function using this regularization 
and the results are almost the same as the ones by the previous strategy. However 
we do not express here. 

Beyond a homogeneous and isotropic model, we next analyze the Bianch type 



I minisuperspace model . The metric of this type can be expressed as 

ds'^ = a'^I^^Lj^dr'^ + a2(r)rfa;2 + b'^{T)dy'^ + c^{T)dz^]. (18) 

The action of this type is 

1 r" [■ a ■ ■ ■ 

I{a,b,c) = - / dr — — {dbc + abc + abc) + NbcX . (19) 

2 J^-f I N J 

In order to deal with the path integral more easily, we change variables of a, b, c as 
X, y, z where x = {bc + a )/2, y = (be — a ), z = a be. The action and the wave 
function are 



/ = / dripxX + pyi) + pzz - NH), 
H = ~pl+pl-^pl-^{x + y), 

^ = dN Dp^DpyDp^DxDyDze^^ = / dN-^{N). 



(20) 



(21) 



The integral with respect to x, y, z is obtained by means of dividing the lapse 
interval N into pieces. Then \E' becomes 



where 



Io = -^{{x"-xf-{y"-yr + 2iz"-zr 



XN 
+ 

Y 



+ ^[{x' + x") + {y' + y")} (23) 

Y 

= XN + 



N 

In this case, the N integration in (22) also diverges, if we take the integration range 
of the A^ from — oo to +oo. Therefore we also use a new regularization method and 
perform the A^ integration. We define again an analytic function by introducing 



a complex parameter a. After the integration, we take the analytic continuation 
q; ^ 1. On the contrary to a homogeneous and isotropic case, we can calculate the 
A^ integration explicitly and the result becomes (up to a normalization factor) 



^ -2-^XY 



v^ = ^-e-^VAr_ (24) 

In Ref. [10], Halliwell and Louko considered the particular Bianch type I model in 
which 6(r) = c(t). Our result may be a part of the results of the contour changing 
method. This reason is that we have treated the integration range of the lapse 
function A^ from — oo to +oo and we have only turned up the straight contour line 
again. By considering a branch cut, we also get a solution that is zero. However 
we do not express above. 

We finally summarize our results. In this paper we have proposed a new regu- 
larization method by concentrating the integration of the lapse function and eval- 
uated the wave function of the Friedmann- Robertson- Walker model along this reg- 
ularization method. Our results are more general because in our method we can 
deal with the branch cut explicitly. We have then evaluated the wave function 
of the Bianch type I model. In this case our result is more limited because we 
treat the integration range of the lapse function from — oo to +oo. While we only 
consider above two simple models, this regularization method will work well in the 
more complicated type or the matter coupled case. 

There remain further works with respect to the interpretation of the wave 
function by this regularization method. Here we compare our method to the con- 
tour changing method. One of the most difficult points in the interpretation of 
the contour changing method is that its form is neither the Euclidean nor the 
Lorentzian path integral ' ' ' . The wave function following this method can't, in 
general, represent the transition from a purely real Euclidean metric to a purely 
real Lorentzian one. This reason is as follows. The contour of this path integral 
is defined as the steepest-descent path with respect to the complex lapse function. 
Thus the contour is a complicated curve. In our method on the other hand, by 



concentrating the lapse function we remain the real Euclidean form and evaluate 
the wave function in a way which is free from a divergence problem. The contour 
of our path integral is defined as the turned up straight line by using the complex 
parameter a in an analytic function. The contour therefore is a simple line. In 
our method, we can calculate the wave function of the universe without using the 
saddle points and the steepest-decsent path. We do not spoil above picture and 
can deal with the real tunneling and get the wave function of the real Euclidean 
path integral more easily. 

The evaluation of the wave function of the universe is very difficult. The 
physical interpretation of the regularization method is not clear. We consider 
there exist some possibilities to obtain a finite wave function along the Euclidean 
path integral regime, and our new regularization method is nothing but one of 
them. 
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